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ABSTRACT 

Observations of the Galactic Centre show evidence of disc-like structures of very young 
stars orbiting the central super-massive black hole within a distance of a few 0.1 pc. 
While it is widely accepted that about half of the stars form a relatively flat disc 
rotating clockwise on the sky, there is a substantial ongoing debate on whether there 
is a second, counter-clockwise disc of stars. 

By means of iV-body simulations using our bhint code, we show that two highly 
inclined stellar discs with the observed properties cannot be recognised as two flat 
circular discs after 5 Myr of mutual interaction. Instead, our calculations predict a 
significant warping of the two discs, which we show to be apparent among the struc- 
tures observed in the Galactic Centre. While the high eccentricities of the observed 
counter-clockwise orbits suggest an eccentric origin of this system, we show the eccen- 
tricity distribution in the inner part of the more massive clockwise disc to be perfectly 
consistent with an initially circular disc in which stellar eccentricities increase due to 
both non-resonant and resonant relaxation. 

We conclude that the relevant question to ask is therefore not whether there are 
two discs of young stars, but whether there were two such discs to begin with. 
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1 INTRODUCTION 

Observations of the Galactic Centre (GC) revealed one 
or two discs of 5-6 Myr old stars orbiting the central 
super-massive b l ack h ole (SMBH) at a distance of ~0.1pc 
l|Pauniard et al.l 120061 . and references therein). The pres- 
ence of such very young stars in the close vicinity of 
the SMBH is puzzling, and a number of proces ses have 
been proposed to resolve this Paradox of Youth: I Gerhard] 
|200ll ) suggested that a disc of stars may have formed 
by tidal disruption of an infalling cluster of young stars, 
which would require a large mass or a central intermediate- 
mass black hole to survive the strong tidal forces from 
the SMBH (iMcMillan fc Portegies ZwartJ [20031 ) . However, 
iLevin fc Beloborodovl (|2003l ) showed that this cannot ex- 
plain the high proximity of the stars, and propose in-situ 
formation by fragmentation of a massive accretion disc a s 
an alternative scena rio (s ee alsolNavakshin fc Cuadrall2005l ). 
iMapelli et ail (|2008l ) and lBonnell fc Ricel (|2008l ) have shown 
the infall of a molecular cloud towards the GC to be effec- 
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tive in creating a disc of stars with distances from Sgr A* 
consistent with the observed young massive stars. 



While the disc of stars orbiting clockwise around the 
SMBH in the centre of our Galaxy is very distinct, the 
disc structure of the counter-clock wise syste m is still un- 
der di scussion l|Lu et al . 2006, 2008). However. lCuadra et al.l 
l|2008h have shown that a single cold disc in the absence of a 
perturbing potential cannot explain the observed large incli- 
nations and eccentricities. The required pertur bing potential 
may result from a massive circum-nuclear disc dGenzel et alj 
1 19941 ). an intermed iate-mass black hole and/or star cluster 
probably IRS 13; iMaillard et al.l l2004h . or a second disc 
Paumard et alj [2006t ). As there is evidence for a second, 
counter-rotating disc, and it has been shown that disc dis- 
solution i s effective enough by m utual interaction between 
two discs (|Navakshin et al.ll2006l ). this is currently the most 
promising solution suggested. We therefore test the hypoth- 
esis that two discs of stars formed almost simultaneously a 
few Myr ago ag ainst the observations (see also review by 
lAlexandeifoorj ). 
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2 NUMERICAL METHOD AND MODELS 

To test the mutual interaction of two stellar discs orbit- 
ing the central SMBH, we perfo rmed direct TV-body sim- 
ulatio ns using our BHINT code (Lo ckmann fc BaumgardU 
2008), which has been developed specifically to calculate the 
dynamics of stars orbiting around a SMBH. It includes post- 
Newtonian approximation of the effects of general relativity 
up to order 2.5. However, test runs show that these correc- 
tions do not have a major impact on the results discussed 
in this work, as the stars hardly reach eccentricities e > 0.9. 
The perio d of relativisti c precession about a 3.5 x 10 6 Mq 
SMBH is l|Einsteinlll9il ) 



fprec.GR = 2 X 10 J P or b 



(l-e 2 ) 



(1) 



v 0.1pc, 

where P or b, a and e are the orbital period, semi-major axis, 
and eccentricity, respectively, while th e period of pre c ession 
due to a sp herical cusp as obse rved bv lSchodel et al.l (|2007l ) 
is given by l|lvanov et al.ll2005l ) 
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Dividing the two yields 

Pprec.GR 



O.lpc 



(1-e 2 )" 



2.7 x 10 J 



O.lpc 



2\15 



(2) 



(3) 



which for a > 0.03 pc (around the minimum achieved for 
the disc stars in our simulations) and e < 0.9 is larger than 
~ 8, showing that precession is dominated by the cusp effect 
in the scope of this Letter (here we will not care about the 
actual distribution of eccentricities e > 0.9), which we also 
find in our simulations. Therefore we waive post-Newtonian 
treatment to keep computing times low. 

Mass ratio and distributions of the stars' orbital param- 
eters are chosen such that the final state is consistent with 
the observed values. We start with two cold stellar discs with 
total masses of 10,000 and 5,000 M Q , respectively, inclined 
at an angle of 130° (run ^4 ). Similar values have been shown 
bv lNavakshin et alj (|2006l ) to be sufficient in leading to an 
extent of disc dissolution comparable to the observations. 
The discs have a finite radial extent of 0.05 — 0.5 pc and 
0.07 — 0.5 pc, respectively, a surface density E(-R) oc R~ 2 5 , 
and a thickness (standard deviation) of 1.44°. The stellar 
orbits have a mean eccentricity of 0.03. The underlying ini- 
tia l mass function (IMF) has a slope of 7 = 1.35, as derived 
bv lPaumard et all (|2006T ) . Since the disc mass for such a flat 
IMF is not very sensitive to the lower stellar mass limit, we 
use a mass range of 1 — 120 Mq . 

The presence of a stella r cusp as it is observed in the 
GC fe.g. Jschodel et alj|2007l) leads to pericentre shift of ec- 
centric orbits, while the orbital plane (and thus any disc 
warping) is conserved. As llvanov et al.l 1 2005h pointed out, 
pericentre shift prevents high eccentricities achieved by the 
Kozai effect. To correctly model the disc evolution, we ac- 
counted for the influence of a background cusp by adding 
14,000 stellar mass black holes (SBHs) of 15 Mq into the 
inner 0.22 pc, consis tent with the density profile derived by 
ISchodel et~atl (120071). The use o f 15 Mq black holes is con- 
sistent with Freitag et alj l|2006l) finding that this region is 
mass-dominated by the most massive remnants having mi- 
grated into the centre by dynamical friction during the life- 
time of the Galaxy. 




Figure 1. Orientation of stellar orbits as a function of semi- 
major axis a after 5 Myr of simulation. The compact shape of 
filled circles indicates moderate warping of the massive disc, while 
the elongated structure of squares shows a massive warping until 
disc dissolution towards the inner edge of the less massive disc. 



Outside 0.22 pc, the eccentricities achieved are only 
moderate, thus the effects of pericentre shift are expected 
to be low. To verify this, we calculated models including an 
analytic potential or a cusp of 40 Mq SBHs outside 0.22 pc, 
which did not yield significantly different results. 



3 RESULTS 

3.1 Disc warping 

The precession frequency of a star orbiting a SMBH of mass 
Msmbh at radius R at an inclination /3 relative to a narrow 
disc of mass Mdia c at a radius Pdisc can be approximated as 
l|Navakshinll2005l ) 



3 M Aisc /GMsmbh 

LOn — r — COS (J' 



4M S iv 



H tt diB 



R 3 



( J R 2 + ^Lc) 5/2: 



(4) 



which for our model roughly translates to a precession pe- 
riod of Pprec.disc ~ 5 Myr (R/0.1 pc)~ 3 ^ 2 for stars in the less 
massive disc (about the massive disc) , and twice as long for 
stars in the massive disc. Note, however, that this is only a 
very rough approximation and is only valid if the disc is not 
radially extended and 7? <C Pdisc- For the evolution of two 
massive discs, the situation is further complicated by preces- 
sion of the reference disc, precession of stars about their disc 
of origin and the resulting increase of |cos/3|, which makes 
detailed predictions impossible. 

As can be seen, the precession frequency is proportional 
to the other disc's mass, but also depends on the distance 
of a star from the central SMBH. Since both discs have a 
finite extent, stars at different central distances precess with 
different frequencies, thus warping the originally flat disc. 

Fig.Q]shows the Aitoff projection of the normal vectors 
to the stars' orbital planes after 5 Myr of simulation. The 
compact shape of the massive disc's normal vectors indi- 
cates the moderate warping of this disc. This also holds for 
the outer part of the less massive disc, while the elongated 
structure at smaller radii shows a massive warping until disc 
dissolution towards the inner edge of this disc. 

The observed GC discs show moderate thickness 
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Figure 2. Orientation of stellar orbits in the GC as a function 
of projected central distance p. Symbols as in Fig. [TJ We find 
a compact spiral structure for the massive (clockwise) disc, and 
an elongated structure depicting the dissolving counter-clockwise 
disc. 



l|Paumard et al.l 120061 ) . To illustrate the present warping, 
we fit discs to subsets of the stars at different radii. Having 
grouped the stars int o a clockwise and a co unter-clockwise 
system according to IPaumard et al.l ((2006), we sort each 
group by projected radius, and fit a disc to the velocity 
vect ors of any 6 subsequent mem bers of the two sorted lists. 

iLevin fc Beloborodovl (|2003h fitted discs to the GC stars 
using a \ 2 test based on the observed velocity vectors and 
three-dimensional error bars. For better comparison between 
observations and simulations, where we have no error bars, 
we fit discs by minimising the root mean square angle a 
between the stars' velocity vectors Vi and the fitted plane, 



1 

N ^ 



arcsin nw. 



(5) 



where n is a normal vector on the fitted plane, a can be 
regarded as the thickness of the disc. To account for obser- 
vational uncertainties, we draw the three components of the 
stars' velocity vectors from a Gaussian using the respective 
observed velocity component and observational error as the 
expectation value and standard deviation, respectively. The 
actual orientation of an observed disc is thus retrieved as 
the mean of the resulting normal vectors from 2000 such 
realisations. 

Fig. [2] shows the fitted normal vectors for the observed 
clockwise and counter-clockwise system. While the massive 
clockwise disc shows a compact spiral structure, demonstrat- 
ing the overall flatness of this disc, the elongated shape of 
the less massive (counter-clockwise) system hints at the fact 
that it loses its disc shape towards the centre, as the orien- 
tation is a strong function of central distance. 

The warp direction in the observations appears to be 
clockwise for decreasing central distance, just opposite to the 
results from our simulation as depicted in Fig. [T] This fact 
suggests that the two discs started with an initial inclination 
of less than 90°. We therefore integrated a second model, 
with an initial inclination of 88° (run B). 

For better comparison with the observations, we now 
sort the massive stars of either disc by semi-major axis, and 
again fit discs to the velocity vectors of now any 12 subse- 



Figure 3. Results from run B with an initial inclination of 88°. 
This time, we plot fitted normal vectors obtained in the same way 
as for the observations. The concordance of direction and shape 
of the elongated structure for the less massive disc is obvious. We 
also find a compact spiral for the clockwise disc similar to the 
observations, although with slightly different orientation. 



quent stars (since we have more stars available than in the 
observations). The results are shown in Fig. [3] Again, we find 
an elongated structure for the counter-clockwise system, this 
time consistent with the orientation of the observed orbits as 
in Fig. [2] The compact image of the clockwise disc exhibits 
a spiral structure similar to the observed disc, although with 
slightly different orientation. The diffusion towards the inner 
edge of this disc is due to fast precession of the correspond- 
ing stellar orbits about their disc of origin. 

To quantify the amount of warping of the stellar discs, 
we can measure the angle between the planes fitted to two 
subsets of stars at different radii within one disc. As we do 
not have the original orientation of the observed discs, we 
take the orientation of the outermost fitted plane as refer- 
ence, since according to Eq. [4] and our simulation results, 
the outermost stars show the least warping. Fig.[4]compares 
the thus calculated warping angles of the observations and 
our run B. It can be seen that the warping is higher for the 
less massive disc, especially towards the inner edge. This also 
explains the lack of observe d counter-clockw ise disc stars to- 
wards the inner edge. Also, iLu et al.l (|2008l ) do not find any 
counter-clockwise disc structure in their improved measure- 
ments within 0.14 pc from Sgr A*, which due to the strong 
warping effect would not be expected, as our simulations 
show. 

Averaging over five realisations of our model, we retrieve 
maximum warping angles of 28° ± 4° for the clockwise and 
108° ± 19° for the counter-clockwise disc. While the warping 
angles of the counter-clockwise disc fit the observations quite 
well, the observed warping of the clockwise disc is higher 
than in our simulations, which may be due to a larger mass 
in the other disc, or any other perturbing potential like star 
clusters or the circum- nuclear disc. 

Similar to what we have seen in Fig. [3] the dips towards 
the inner edge of the massive clockwise discs can be ascribed 
to fast precession of the corresponding orbits about their disc 
of origin. 
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Figure 4. Disc warping angle as a function of central distance. 
We plot the angle between the plane fitted to any disc's subset of 
stars and the outermost fitted plane for the respective disc. The 
upper panel shows the results for the clockwise disc from both ob- 
servations (red) and simulations (run B\ green); the lower panel 
shows the respective results for the counter-clockwise systems. 
Horizontal and vertical error-bars indicate standard deviation of 
central distance and lcr-thickness of the fitted planes, respectively. 
It is evident that the angle does not randomly scatter, but in- 
creases with decreasing central distance, thus showing the disc 
warping. 




3.2 Relaxation in the Galactic Centre 

As we have shown in the previous section, the precession 
caused by the stellar discs is the main driver for the change 
of the stars' orbital inclination. However, it does not change 
orbital eccentricity or semi-major axis. These properties are 
changed by relaxation processes among the stars in the GC. 

The relaxation time of a stellar system is known as the 
time it takes for the stars to change their energy (or angular 
moment um) by a fact or of order unity and can be estimated 
as (e.g., ISpitzerlll987l . Eq. (2-62)) 

_ 0.065a 3 

— G 2 M*p In A' W 

where a is the local 3D velocity dispersion, p is the stellar 
mass density, and In A is the Coulomb logarithm. Using the 
orbital circular velocity aro und a 3.5 x 10 6 Mq SMBH for 
a, the stellar mass profile bv lSchodel et al.l ((2007) at central 
distance r < 0.22 pc, and In A ~ 10, this gives 



t le i ~ 190 Myr 



0.1 pc 



M, 

15 Mr. 



(7) 



The mean squared change of energy E and angular momen- 
tum L of a star after a time t is described by 

t/t le i = ((AE/Eof) = ((AL/L c ,o) 2 ) (8) 

ijRauch fc Tremainelll996l eq. (4)), where Eq is the initial 
energy and L Ci o the angular momentum of a corresponding 
circular orbit. This way, we can expect the stellar orbits in 
the GC to be fully randomised only after ~ 200 Myr, and 
only some stars on moderate eccentricities after 5 Myr. How- 
ever, relaxation is enhanced within the stellar discs: While 
the stellar discs do not contribute much to the average den- 
sity of the cusp, they significantly increase the local density 



in the vicinity of the disc stars. In our models, a sphere of 
radius r 3 <C 0.1 pc centred on a star in the massive disc 
at 0.1 pc from the SMBH contains 3 x 10 7 (r 3 /pc) 3 M of 
cusp stars and 2.6 x 10 5 (r s /pc) 2 Mq of disc stars, hence for 
r s ~ 0.01 pc the disc stars double the local density and thus 
the amount of relaxation. The small relative velocities in the 
initially cold disc are expected to further enhance this effect. 

Fig. [5] shows the evolution of mean squared change 
of energy and angular momentum for the stars of run B 
for different initial semi-major axes. While the cusp stars 
follow the theoretical prediction of energy relaxation very 
well, it is clearly seen that relaxation is faster for the disc 
stars, especially so in the beginning, where the local den- 
sity is at its maximum. Furthermore, angular momentum 
relaxation appears to be much faster than predicted by 
non-resonant relaxation and shows a strong dependence on 
central distance, suggest ing that resonant relaxation (RR; 

iRauch fc Tremainell 19961 ) plays an impo rtant role. 

The RR time-scale is given by IRauch fc Tremainel 

1 19961 . eq. (9)) as 



-MsMBH 



370 Myr 



0.1 pc 



15 Mr 



(9) 

hence, while the contribution of RR at 0.1 pc is only half as 
much as for non-resonant relaxation (NR), it is expected to 
dominate relaxation further in. Replacing the stellar mass 
of M* = 15 Mq by the Paumard mass function used for the 
discs would increase RR efficiency by a factor of 3, since M, 
needs to be replaced by J m 2 dN{m)/ J m&N(m) ~ 50 Mq. 
The alignment of stellar orbits in two discs may further ac- 
celerate RR, as the sum of torques exerted causing a change 
in angular momentum is no longer random. Replacing the 
stellar cusp by an analytic potential, we find that angular 
momentum relaxation as shown in Fig. [5] is decreased by 
~ 20%, but still significantly higher than that of the cusp 
stars in the previous simulation, confirming that while the 
cusp stars do contribute, relaxation among disc stars is dom- 
inant. 

Note that RR is only fully effective as long as the orbits 
are almost stationary over some time t such that P or b <C t <C 
fprcc, which according to Eq. [5] is only fulfilled in the inner 
part of the discs; however, in the outer parts of the disc the 
theoretical RR time-scale is much longer than the NR one, 
so that resonant effects cannot be expected to significantly 
contribute to stellar relaxation. 

In the following, we analyse the distribution of the disc 
stars' orbital parameters resulting from relaxation as de- 
scribed above. Fig.[S]shows the eccentricities of massive disc 
stars (> 17 Mq) after 5 Myr of simulation as a function of 
semi-major axis. It can be seen that a considerable fraction 
of stars from both discs achieve eccentricities up to e = 0.9. 
This shows that the observation of young stars on highly ec- 
centric orbits is as such not a valid argument against a circu- 
lar disc origin of these stars. However, this is only true for the 
inner part of the disc, while most stars beyond 0.3 pc from 
the SMBH do not show eccentricities larger than e = 0.2 
in our runs. Since this might be explained by the absence 
of cusp stars in this region, we computed another model in- 
cluding a stellar cusp extending beyond the outer edge of the 
discs, which did not yield significantly higher eccentricities. 

Fig. [7]shows the eccentricity distribution of stars within 
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Figure 5. Time evolution of mean squared energy and angu- 
lar momentum change for model B, averaged over five different 
realisations with standard deviation shown as error bars. The 
normalisation constants Eq and L c ,o are the local energy and an- 
gular momentum of a circular orbit, respectively. It can be seen 
that relaxation in the discs (dotted lines) is well enhanced, most 
likely due to local over-density and small relative velocities. En- 
ergy relaxation (left) of the cusp stars follows the theoretical pre- 
diction of non-resonant relaxation very well, showing only weak 
dependence on central distance. The apparent radial dependence 
of energy relaxation in the discs may be explained by faster de- 
crease of local over-density as a consequence of stronger warping 
at the inner edge of the disc (dotted red line). On the other hand, 
relaxation of angular momentum (right) shows a strong radial 
dependence and appears to be much faster than energy relax- 
ation towards the centre. We thus conclude that an additional 
component of angular momentum relaxation must be present, i.e. 
resonant relaxation. 
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Figure 6. Eccentricity of massive disc stars as a function of (pro- 
jected) central distance. The filled circles show the results of run B 
after a simulated ti me of 5 Myr, the open sq uares depict the obser- 
vations as given by iPaumard et al.l [|2006f ) . Observational errors 
are denoted as la error bars. While the eccentricities observed 
in the clockwise disc (left panel) within a distance of < 0.3 pc 
from the SMBH can be well explained by our model, the orbits of 
stars towards the outer edge of the population (having an admit- 
tedly high uncertainty) may need a different explanation, such as 
a more complex infall geometry, or a massive perturber. 
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Figure 7. Eccentricity distribution of stars within 0.3 pc from 
the SMBH. The green and cyan boxes show the results of our 
models A and B, respectively, each averaged over five different 
realisations with standard deviation sh own as error bar s . The red 
filled boxes depict the values derived bv lPaumard et all d2006T l for 
the observed discs with Poisson error bars, (a) The eccentricity 
distribution of the observed clockwise disc is in good agreement 
with our simulations, (b) In the observed counter-clockwise sys- 
tem, 6 out of 11 stars have eccentricities e > 0.7, which is not 
consistent with our simulations. 



0.3 pc from the SMBH. While the observed eccentricities of 
the clockwise disc are very well resembled by our simula- 
tions (upper panel), the fraction of 6 out of 11 stars with 
observed eccentricities e > 0.7 in the counter-clockwise sys- 
tem is significant and not consistent with any of our two 
models. However, integrating a model where the less mas- 
sive disc starts with the observed eccentricity distribution 
results in a consistent eccentricity distribution of this disc 
while retaining its warping and the other disc's properties. 
We thus conclude that the stars in the counter-clockwise sys- 
tem prob ably formed from an ec centric disc, which has been 
shown bv lBonnell fc Ricd (|2008l ) to be a viable scenario. 



4 DISCUSSION AND OUTLOOK 

Our results show that the orbits of young stars observed 
within 0.3 pc from Sgr A* can be very well explained by 
the interaction of two initially flat discs. While the present 
existence of a counter-clockwise disc in the GC is still con- 
troversial, the observed structures do not exclude but natu- 
rally result from a setup of two inclined accretion discs. We 
find that this scenario especially accounts for the controver- 
sial nature of the counter-clockwise disc (i.e., a higher disc 
thickness) , as this less massive system is exposed to stronger 
warping and dissolution. Especially towards the inner 0.1 pc, 
this warping is so strong that one cannot expect to find a 
disc structure after 5 Myr of i nteraction. Hen ce, the two- 
disc scenario is consistent with iLu et al.l l|2008l ) not finding 
a counter-clockwise disc in their recent data set extending 
to 0.1 pc. Whether or not the less massive stellar system is 
called a disc may be a matter of discussion, even if it initially 
had a perfect disc shape. 

From our simulations we see that both non-resonant 
as well as (short-term) resonant relaxation effects play a 
significant role in the evolution of orbital eccentricities of 
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the disc stars, although the exact contributions are not yet 
fully understood. Furthermore, we see that both disc and 
cusp stars contribute significantly to relaxation. 

We are left with the question of how two mutually 
highly inclined dense accretion discs can have formed (al- 
most) simultaneously a few Myr ago. Furthermore, our mod- 
els do not explain the relatively large eccentricities and in- 
clinations of the orbits towards the outer edge of the stellar 
disc(s). For this, one may consider the effects of an external 
perturbation, or simply argue that the uncertainties of the 
obse rvations of these stars are comparably high. In addi- 
tion, lBonnell"fc Rice! | 20081 ) have shown how eccentric stel- 
lar discs can be created from high-eccentricity infalls of giant 
mol ecular clouds (GMCs ) . 

iHobbs fe Navakshinl (|2008l ) simultaneously explained 
both the double accretion disc and the diffusion at larger 
radii: Their calculations show how (self-) colliding GMCs get 
bound to the SMBH and may form two accretion discs at 
large inclination in the vicinity of the SMBH. At somewhat 
larger distances, the geometry is much more complex, which 
may account for the observed larger eccentricities and incli- 
nations. 

We suspect that a thorough analysis of the orbital pa- 
rameters of young stars in the GC will contribute signifi- 
cantly to our understanding of star formation under such 
extreme circumstances. Our calculations rule out the exis- 
tence of two 5 Myr old flat discs within the inner 0.3 pc from 
Sgr A*. None the less, detailed properties of the stellar or- 
bits as obtained e.g. by analyses of stellar accelerations may 
help to describe the disc warping and eccentricity distribu- 
tion more precisely and thus shed light on the exact gas 
infall scenario. Furthermore, the detection of B-type stars 
in the discs is expected to increase the significance of the 
disc warping (and of the counter-clockwise disc, as larger 
data sets may yield the opportunity to confirm it despite 
massive warping by constraining to subsets of smaller radial 
extent) and yield more information about the discs' stellar 
mass function and total mass. 
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